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Abstract 



In a quantum system with a smoothly and slowly varying Hamiltonian, which 
approaches a constant operator at times t ±00, the transition probabilities 
between adiabatic states are exponentially small. They are characterized by 
an exponent that depends on a phase integral along a path around a set of 
branch points connecting the energy level surfaces in complex time. Only cer- 
tain sequences of branch points contribute. We propose that these sequences 
are determined by a topological rule involving the Stokes lines attached to 
the branch points. Our hypothesis is supported by theoretical arguments and 
results of numerical experiments. 
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I. INTRODUCTION 

Consider a quantum mechanical system with a Hamiltonian H{X) that depends analytically 
on a parameter X and has a discrete spectrum. Suppose the parameter changes slowly and 
analytically with time t, in such a way that asymptotically as t — ^> ±oo, it remains constant. 
It is convenient to write X = X{t) where r = et. The adiabatic parameter e is a small, real 
number which provides a dimensionless measure of how slowly the Hamiltonian changes, 
due to the change in the parameter X. The quantum adiabatic theorem states that if such a 
system is prepared in the rfi^ eigenstate of the instantaneous Hamiltonian at any time, then 
in the limit e — 0, the system remains in the n**^ eigenstate of the instantaneous Hamiltonian 
at later times [1,2]. 

Corrections to the adiabatic theorem can be extremely small. If X{t) approaches its limiting 
values sufficiently rapidly as r ±oo, the probability of making a transition to another 
eigenstate of the Hamiltonian vanishes exponentially (faster than any power of e) as e — *• 0. 
The probability for making a transition from the state |0„(X(ro))) to the state \(l)m{X{Tf))) 
is given asymptotically in e by 

Pn^m ~ C„„exp[-25'„„//i] (1.1) 
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where Cnm and Snm are positive constants and the action Snm is inversely proportional to 
e. Results of this type have been proven rigorously for most cases of physical interest in 
which the Hamiltonian is a 2 x 2 [3,4] or 3 x 3 matrix [5,6], and for some special cases of 
multi-level systems [7-10]. This paper will discuss how these constants Snm, describing 
non-adiabatic transitions may be determined with some generality for a discrete multi-level 
system. For definiteness, we will assume that the Hamiltonian if is a symmetric J\f x Af 
matrix, which is real for real r, and analytic on a strip containing the real r axis. Then 
the Schwartz reflection principle guarantees conjugacy of the matrix elements about the real 
time axis. 

The special case of a two-level system has been discussed by many authors. Zener [3] 
calculated an exact expression for the transition probability for a two level model of the 
form 



This Landau-Zener system is of particular importance because for a multi-level system, 
according to degenerate perturbation theory, a near- degeneracy or 'avoided crossing' of two 
levels can be approximated by a two-state Hamiltonian of this form. The exact transition 
probability for this system is 



Note that the closest approach of the energy levels E{t) is A, so that non-adiabatic transi- 
tions occur easily if the energy levels approach each other closely. 

Later Dykhne [4] showed that in the hmit e — > the transition probability for a general 
two level system is determined entirely by the analytic continuation of the energy levels as 
a function of complex time. The eigenvalues of the instantaneous Hamiltonian E{t) form 
two branches of a single complex function, which are connected at branch points occurring 
as complex conjugate pairs. The action exponent S is determined by integrating the energy 
along a path which leaves the real axis and loops around a branch point so that it returns 
on the other branch. Dykhne showed that, in the limit e — > 0, the transition probability 
becomes 



where 7 is the path illustrated in figure 1. The pref actor in (1.1) is seen to be C = 1. If 

there arc more than one pair of branch points, the one with the smallest value of S gives the 
dominant contribution to the transition probability. For the model (1.2), S — nA'^/AeA, so 
that the Landau-Zener formula follows from this more general result. 

In a multi-level system, it is natural to expect that the transition probability is, in general, 
the product of several terms of the form of (1.4), corresponding to several successive transi- 
tions between pairs of levels. The conclusions of this paper will confirm this expectation. It 
might also be expected that the transition probability will be determined by the combination 
of paths around branch points that minimizes the total value of the sum of integrals S of 
the form appearing in (1.4). This expectation is only partially correct: not all sequences of 
branch points are allowed. In this paper we propose a topological constraint, or rule, that 
the allowed combinations of branch points should satisfy. The value of S is determined by 




(1.2) 



Plz = exp(-7rA72Ae;i) . 



(1.3) 




(1.4) 
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the minimal value of the sum of integrals of the form (1.4) for paths consistent with the 

constraint. 

An heuristic reason for the necessity of a rule that eliminates certain sequences of branch 
points can be understood from an example given by Hwang and Pechukas [5]. Consider a 
three level system, with energy levels as plotted schematically in figure 2(a). There are two 
avoided crossings where pairs of levels approach each other very closely. In the neighborhood 
of these avoided crossings the two nearly degenerate levels can be represented by a two-state 
model similar to (1.2), and the transition probabilities between pairs of adjacent levels are 
determined by the branch points near the real axis shown in figure 2(b). The transitions 
from level 1 to level 3 occur through Landau-Zener transitions at the successive avoided 
crossings, and the transition probability is the product of two terms of the form of (1.4). 
The transition from level 3 to level 1 cannot make use of the same path, however, because 
the avoided crossing between levels 2 and 1 is encountered before the transition from level 
3 to level 2 has taken place. Instead, the transition occurs through a more distant branch 
point between levels 1 and 3 as shown in figure 2(b). The probability for the 3 — > 1 transition 
is therefore expected to be much smaller. Therefore, a rule is required that eliminates the 
path from level 3 to level 1 via the branch points closest to the real axis, but permits this 
path for the reversed transition. 

In section 2 a rule is proposed for determining which branch points contribute to the transi- 
tion probabilities. Our rule is motivated by a discussion of the Stokes phenomenon given by 
Heading [11,12], which assumes that evanescent waves are meaningful even in the presence 
of larger contributions, and that the amplitudes of evanescent waves change on the Stokes 
lines. These assumptions must be justified by a reference to an asymptotic series in e for the 
wavefunction, showing that the errors in the dominant terms can be reduced to below the 
magnitude of the subdominant (evanescent) terms, except in the vicinity of the Stokes lines. 
This interpretation of the asymptotic series for the wavefunction occurs in Stokes' analysis 
of the Airy function [13], and many further examples are discussed in a book by Dingle [14]. 
Works by Berry [15-17] have developed the theory further and discuss physical applications, 
including quantum adiabatic theory for two-level systems. 

The standard approaches for this interpretation of Stokes fines are vafid only for those Stokes 
lines associated with branch points having the smallest action S, since these give the most 
divergent contributions to the asymptotic series. All applications to date fiave considered 
systems with only two states (in adiabatic theory) or two channels (in WKB approaches 
to scattering problems), where only this one branch point is relevant. These methods are 
bound to fail in multi-level systems where many branch points are of necessity involved in 
a transition, since the divergence of the series is dominated by the branch point with the 
smallest action. In section 3 we describe a formulation that allows all of the Stokes lines to 
be given a similar interpretation. By projecting out tfie subspace corresponding to brancfi 
points witfi small actions, divergences due to otfier brancfies are revealed. Tfie metfiod is 
based upon an exact iterative 'renormalization' of tfie Hamiltonian, adapted from a scfieme 
introduced by Berry [18]. Section 4 discusses tfie limiting form of tfiis series of renormalized 
Hamiltonian operators and its interpretation in terms of tfie Stokes pfienomenon. Section 
5 contains some numerical illustrations of tfie results, and section 6 summarizes tfie new 
contributions and considers furtfier development of tfie tfieory. 
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II. A TOPOLOGICAL RULE FOR DETERMINING THE ACTIONS 



In this section we propose a topological rule for determining the actions for non-adiabatic 
transitions motivated by a discussion of the Stokes phenomenon. The first three subsections 
review some ideas well known in the context of WKB theory. We review these ideas here to 
place them in the slightly different context of adiabatic theory and to establish our notation. 

A. Adiabatic solutions of the Schrodinger equation 

It is possible to write down an asymptotic series for the solution of the Schrodinger equation 
ihdt\^p{t)) = H{t)\ilj{t)) as an expansion in powers of e 

\Mt)) = eM-iOn(t, io) A]{|0n(T)> + e|(/>«(T)) + ... + e'^l^f (t)) + ...| (2.1) 

where |</'ra(T)) is an eigenvector of the instantaneous Hamiltonian H{X{t)), and 9n{t,to) is 
the phase integral of the corresponding eigenvalue (n*^ energy level) defined by 

Onit, to) = f dt' E^{t') = - r dr Enir) , (2.2) 

J to ^ Jto 

with ^0 the phase reference. The state in (2.1) is an approximate solution of the time- 
dependent Schrodinger equation, assuming the system was prepared in the n*'^ energy level 
in the remote past, that doesn't take into account non-adiabatic transitions. Only the 
leading order terms in these asymptotic expansions will be considered explicitly. These will 
be termed the adiabatic solutions 

\Mt)) = eM-^On{t,to)/h]\Mt)) . (2.3) 

For complex t, the exponential factor need not have unit modulus. As we move along a 
path in the complex t plane, we can compare the rate at which two solutions increase or 
decrease exponentially. The solution for which the derivative of Im 6n is largest is said to be 
dominant and the other solution suhdominant. 

Because the solutions are to be considered as functions of a complex time variable, it is nec- 
essary to include a few comments about how the eigenvalues and eigenvectors are continued 
into the complex plane. Since for real t the matrix elements of the Hamiltonian are real, we 
define the inner product of two eigenvectors without complex conjugation. We shall assume 
that the eigenvectors are normalized in the usual way, using this inner product. As t be- 
comes complex, the eigenvalues and eigenvectors also become complex, but the Hamiltonian 
matrix remains symmetric and its eigenvectors orthogonal. By choosing the eigenvectors 
real on the real axis, they then have no gauge freedom. They may be expressed as analytic 
functions of t at all points in the complex plane, except for a set of singular points. 
Singular points of the eigenvalues and eigenvectors can arise either because the matrix ele- 
ments themselves have singularities, or because of singularities in the mapping from matrix 
elements to eigenvalues and eigenvectors. This latter type of singularity is associated with 
points where eigenvalues become degenerate in the complex t plane. These singularities 
arising from degeneracies will be termed branch points, because they connect two branches 
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of the Riemann surface of the eigenvalue function E{t). The former type of singularities are 
of less interest, because they are not universal in form, and because they typically lie farther 
into the complex plane than the branch points. (The reason for this will be discussed in 
section 5.) Degenerate perturbation theory shows that in the neighborhood of two eigen- 
values becoming degenerate, the energy difference is the square root of the discriminant of 
the corresponding quadratic equation. At the singularity, the discriminant vanishes while, 
generically, its time rate of change is non-zero. Thus generically, singularities due to eigen- 
value degeneracies are square-root branch points. In this paper we shall assume all relevant 
singularites are of this form. 

For systems with a single parameter X, assumed here, a theorem of Von Neumann and 
Wigner [19] states that, generically, degeneracies do not occur for real t. Therefore the 
eigenvalues and eigenvectors can be labeled with an index corresponding to the ordering of 
the (real) eigenvalues on the real axis. Usually it will be useful to consider the eigenvalues 
and eigenvectors to be single- valued functions of t. This requires the introduction of branch 
cuts that restrict the domain of definition of the function E{t) to a single Riemann sheet. 
Typically we will choose the branch cuts to be lines with constant Ret that do not cross the 
real illustrated in figure 2(b). 

It is necessary to establish how the eigenvalues and eigenvectors change upon crossing a 
branch cut. Consider a branch point due to the degeneracy = En2- Since this is a square- 
root branch point, the eigenvalues are simply exchanged as the branch cut is traversed. The 
eigenvectors are more involved. If, starting from a point t+ adjacent to the branch cut, 
the eigenvector |0„^(t+)) is followed anticlockwise around the branch point, upon reaching 
the other side of the branch point the state |0„j(t_)) is equal to a multiple of the state 
|0„2(^+))- This is illustrated in figure 3(a). Upon taking this multiple of the state |0„2(t+)) 
anticlockwise to it is transformed into — |0„^(t_|_)), i.e., the eigenvector recurs with its 
sign changed after two circuits around the degeneracy. This fact is easily verified in the 
special case of the Landau-Zener model, (1.2). Choose a branch cut crossing the real axis 
connecting the complex conjugate branch points, and consider a circuit constructed from the 
real axis and from a semicircle at infinity, as illustrated in figure 3(b). Because the circuit 
can be shrunk to a circle enclosing the branch point, and because the singularities of this 
model are generic, the result is true in the general case. Thus if we choose the phase of the 
state 1 0^2) relative to \(pni) appropriately, we may write 

|0„,(t_))=i|0„2(^+)) 

|0„2(M) = i|0n,(i+)) (2.4) 

as a generic description of how the eigenvectors change when crossing a branch cut in the 
anticlockwise sense. 

B. Solutions in the neighborhood of a branch point 

This section considers the solution of the Schrodinger equation in the vicinity of a branch 
point. The arguments are an adaptation of those given by Heading [11,12] for the case of 
WKB approximations. They arc included here because the terminology must be re-defined 
for the adiabatic problem, and because there is a minor difference in the logic. 
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We assume that at f* the levels E^^ and become degenerate. In the neighborhood of 
t* the difference between the two nearly degenerate levels behaves as AE = En^ — E^^ 
const. (t — t*)^/^, so that the phase integral difference has a singularity of the form 

A0{t, t*) = /* dt' AE{t') ^^{t- t*)=^/2 (2.5) 

Jt* e 

where Kt* is a complex constant characterizing the branch point at t* . Thus it is also 
necessary to introduce branch cuts in the phase integrals 9n{t,t*), which we take to be the 
same as for the energies. On crossing a branch cut we must change the labeling of the energy 
levels and phase integrals, although their values change smoothly. 

Some of the level curves (contours) of the phase integral that pass through the branch point 
are very important in the formulation of the theory. These are the Stokes and anti-Stokes 
lines defined by 

Re A6I = Stokes lines 

Im = anti — Stokes lines . (2-6) 

On the anti-Stokes lines the solutions (2.3) connected by the branch point are co-dominant, 
and on the Stokes lines one solution can be regarded as being maximally dominant over the 
other. The Stokes lines associated with a typical branch point are sketched in figure 4. 
The subdominant solution is meaningful only if it is larger than the error term of the 
dominant series, which can itself be assumed to be comparable to the smallest term of the 
asymptotic expansion (2.1). On considering how the solution behaves near a branch point, 
it is convenient to first consider the behavior of the solution on the anti-Stokes lines. This 
is because the adiabatic series corresponding to the two levels which become degenerate at 
the branch point are co-dominant on the anti-Stokes lines, and are therefore certainly both 
meaningful there. If these solutions are |^i(t)), |^2(^)), fbe solution at any point along the 
anti-Stokes line (Al) in figure 4 may be written 

m))-<^i\Mt))+(y2\Mt)) (2.7) 

with the multipliers ai, a2 approximately constant. On the other two anti-Stokes lines the 
solution takes the same form, but it need not have the same coefficients; for example, on 
line (A2) the coefficients could be [a[,a2). Since the Schrodinger equation is linear, there 
must be a linear relationship between the coefficients on the anti-Stokes lines: 

«n=Mf"0- (2-8) 

If iV'i(^)) is the dominant solution in the sector between the lines (Al) and (A2), then 
the multiplier of this solution cannot change, i.e., a[ = a\. Heading [11,12] has given a 
very general argument, showing that the multiplier of the subdominant solution is altered 
on crossing the Stokes line. For the generic case, where the singularity is given by (2.5) 
(triplets of Stokes lines attached to a branch point) , the argument follows closely that given 
by Heading for the WKB approximation and the result is the same: the transition matrix 
from the anti-Stokes hue (Al) to (A2) is 

^=(-i l)- ^2-^) 



7 



In other words, — i times the multipher of the dominant solution is added to that of the 
subdominant solution whenever a Stokes line is crossed in the anticlockwise sense. This is 
Heading's rule. It assumes a phase reference at the branch point. Similar relations exist on 
crossing the other Stokes lines. However, in writing down expressions for these relations it 
should be remembered that, on crossing anti-Stokes lines, it is possible for the dominance 
of the solutions to switch. For example, if \'ip2{t)) is dominant on both (S2) and (S3), 
then dominance switches upon crossing (Al) and (A2) but not (A3). The exact solution of 
the time-dependent Schrodinger equation must be analytic and single valued at the branch 
point, so equation (2.9) is obtained by requiring that the multipliers ccj return to their 
original values when traced in a circuit around the branch point, crossing three Stokes lines 
and the branch cut. The branch cut is accounted for by noting that, according to (2.4), 
upon crossing the cut the labels of the states are exchanged and the states are multiplied 
by i (this circumstance differs from the case of WKB theory, where the factors of i have a 
different origin). We can introduce the matrix 



to account for the switching of dominance on the anti-Stokes lines. The branch cut is 
accounted for by the matrix iT. The effect of making a circuit around the branch point 
is then described by the product iTMTMTM, which is the identity matrix, verifying that 
(2.9) is consistent with a single- valued wavef unction. 

It is desirable to comment on the interpretation of (2.9). The adiabatic wavefunctions are a 
poor approximation to the exact solution at the branch point because they are singular there. 
They are also a poor approximation far away from the branch point because of the presence 
of other branch points and singularities. Furthermore, their interpretation is ambiguous at 
all points off the anti-Stokes lines, because the subdominant wavefunction may be smaller 
than the error of the dominant wavefunction. As they stand then, the arguments above 
suffice only to discuss the behavior of the wavefunction in a (deleted) neighborhood of the 
branch point, near the anti-Stokes lines. 



The results of section 2.2 are sufficient to enable the transition probability due to a sin- 
gle branch point to be determined, provided that they are supplemented by an additional 
assumption. The logic of the argument presented in section 2.2 gives the form of the wave- 
function on two of the anti-Stokes lines, provided it is known on the third (and even on 
the anti-Stokes lines it is only a good approximation sufficiently far from the branch point). 
The additional assumption is that the solutions may be extended from the anti-Stokes lines 
as required and still remain meaningful. This assumption requires the use of asymptotic 
series approximations, which reduce the errors of the dominant terms until the subdominant 
terms are meaningful. Such approximations are not discussed explicitly in Heading's work 
and must be verified for specific appHcations. 

It will be convenient to introduce the following notation for the phase integral factor in the 




(2.10) 



C. Transition probabilities deduced from Heading's rule 



8 



adiabatic solutions: 



fn{t,to) = exp[-i9n{t,to)/h] = exp 



^ fdt'E^it') 
n J to 



(2.11) 



We assume that the system begins in the state \i}i{t)) = to)|0i(^)) as ^ — > — oo. The 
phase reference is an arbitrary point on the real axis. As t — > oo, the system evolves to the 
state |'0i(^)) + 0'\i-'2{^^): where a is the transition amplitude that we seek. In the domain in 
figure 4 bounded between (Al), (SI) and the real axis, the system is in the dominant state 
with respect to the branch point in the upper-half plane. Heading's rule, given by (2.8) and 
(2.9), gives the evolution of the solution around the branch point from (Al) to (A2) as 

/i(i,r)l0i(i))^/i(i,r)l0i(i))-i/2(i,ni</'2(^)) , (2.12) 

or in terms of the adiabatic solutions with phase reference as 

/i(to,r)|V^i(^)) ^ /i(to,r)|^i(t)) -i/2(to,t*)|V^2(t)) . (2.13) 

Assuming we may extend this solution down to time t on the real axis, the probability for 
the system to be found in state |'02(^)) is 

Pi^2 ~ |/2(to,t*)/i(t*,to)|' = exp[-21m5i,2/;i] (2.14) 

where 

^1,2 = /* dt' [E2{t') - Ei{t')] = - /' dt'E{t') (2.15) 

J to 

is the action evaluated around a path, analogous to that shown in figure 1, encircling the 
branch point at t*. Its value is independent of the phase reference to. Thus Dykhne's formula 
(1.4) has been deduced from Heading's rule. 

Note that the branch point below the real axis does not make any contribution to this 
transition probability, because the solution lipi) is subdominant in the lower half plane. To 
compute the transition probability for transitions in the opposite direction, the roles of the 
two conjugate branch points are reversed. 



D. The case of two branch points 

In this section we discuss the case in which a transition occurs in two stages, first from the 
level no to level rii, and then from level rii to level n2. Consideration of this case allows us 
to identify a rule for determining the circumstances under which successive transitions may 
occur. 

We examine a case in which adjacent pairs of levels become degenerate at distinct branch 
points. If we are concerned with transitions from lower to higher energies, branch points in 
the upper half plane are relevant, because for paths below these branch points lower lying 
levels are dominant. The Stokes lines (defined by (2.6)) of the two branch points can assume 
two possible arrangements as illustrated schematically in figures 5(a),(b). 
First consider the arrangement in figure 5(a), with the initial state being such that only the 
level with index no is occupied ast — > — oo. Following the argument used in the last section. 
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we use Heading's rule to connect the solutions on the anti-Stokcs lines in the vicinity of 
the branch points, and assume that these solutions may be extended down to the real time 
axis. Since only the initial state |V'no(^)) occupies the sector between the lines (Al) and (SI) 
attached to the branch point (no,ni), and since lipnoi't)) is the dominant solution on (SI), 
Heading's rule gives the wavefunction on (A2) as a multiple of \ipno{t)) Thus the 

evolution of the solution around the (no,ni) branch point from (Al) to (A2) is given by 

/no(^,C,nJl0no(t)) - /no (t, ^^o.n J I'/'no (t)) " i/n, (t, C„ J |0n, (t)) • (2.16) 

Next, since no Stokes lines are intervening, we may extend this solution upward to the 
anti-Stokes line (Al') connected to the (ni,n2) branch point. The sector between (Al') and 
(SI') is now occupied by the solution \ipni{t)). Since this is the dominant solution on (SI'), 
Heading's rule gives the wavefunction on (A2') as a multiple of \tjjni{t)) — i|^„2(^))- The 
evolution of the solution around the (711,^2) branch point from (Al') to (A2') is given by 

Combining these two evolutions, and expressing the result in terms of the adiabatic solutions 
with a phase reference to on the real axis, one obtains 

l^n„(t)) ^ IV'noW) -i/ni(^0,Co,ni)/no(Co,ni,^o)|^ni(t)) 

U)\^l^nM (2.18) 

for the evolution of the wavefunction from (Al) to (A2'). Extending this solution down to 
time t on the real axis, one finds the probability for making the transition from level Uq to 
level n2 is 

Pno^n2 ~ l/n2(^0? ^rii,n2)/ni(^ni,n2' ^no,ni)/"o(^no,ni' ^o)| 
m.n-z' ^o)/rii(^0) ^no,ni)/no(^no,ni' ^o) | 

= exp[-2 Im {Sn,,m + Sn,,n,)/h] (2.19) 

where Sno,ni and S'„j,n2 are actions for paths evaluated around the branch points located at 
times f^g^m and tj^^^^j, respectively. The transition probability is therefore the product of 
two factors of the Dykhne formula form, corresponding to two successive transitions. 
Now consider the arrangement in figure 5(b) with the same initial state \ipno{t)) occupied as 
t — > — oo. At the (711,712) branch point, |'0„^(t)) is dominant and {ipnii't)) is subdominant. 
According to Heading's rule, a transition into the n2 level requires that the subdominant 
solution be switched on in proportion to the dominant solution's multiplier. However, the 
solution iV'ni(^)) has multiplier zero in the entire region containing the relevant anti-Stokes 
lines attached to the (ni,n2) branch point. It follows that no transition from no to 712 is 
possible using these two branch points. 

There usually will be other branch points farther out in the complex plane connecting levels 
no and n2 directly. In both of the cases discussed above, this branch point would make a con- 
tribution to the transition probability of the form P„g_>„2 = -fn2->no — 6xp[— 2 Im5'no,r^2/^]• 
In the case of figure 5(b), this would be the only contribution. In the case of figure 5(a), it 
is the dominant contribution if lmSno,n2 < -|- lmSni,n2J and negligible otherwise. 
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E. A criterion for selecting possible transition sequences 



We now propose a rule for determining the combinations of actions that correspond to 
allowed transitions. The rule is based upon Heading's local analysis of the form of the 
solution in the vicinity of a branch point, supplemented by the assumption that solutions 
may be continued away from the branch points. 

In the case of an upward transition from level no to n^, a transition might be possible using 
a sequence of branch points Ci,n2> Kik-i,nk with transition probability Pn^^nk — 

exp[— 2Im (5'„o,„i + .. + Sni^_^,nk) /^]- This transition sequence is allowed only if the branch 
points and their associated Stokes lines satisfy a topological criterion. The transition Uj — > 
rij+i is mediated by the branch point at , and can occur only if the level Uj is occupied 

at this branch point. This requires that the branch point t* . ^^ .^^ lies in the quarter plane 
above the real axis, and to the right of the boundary formed by the Stokes line descending 
from t*^._^ ,^. to cross the real axis, and the branch cut from the t*^._^ ,^. branch point. 
The rule is topological in character. For downward transitions, the relevant branch points 
are in the lower-half plane, but because of reflection symmetry, the rule can be applied in 
exactly the same way as for upward transitions. The rule might at first sight appear to have 
a degree of arbitrariness, in that it refers to the positions of the branch cuts as well as the 
Stokes lines. It may, however, be verified that moving the branch cuts does not affect the 
predictions. The reason is that when a branch cut from the t*^^^^ branch point is moved 
past the „2 branch point, the labeling of the levels must be changed so that the latter 
branch point now connects levels Uq and 712 rather than ui and 712. 



III. ADIABATIC RENORMALIZATION WITH PROJECTIONS 

A. Motivation for renormalizing the Hamiltonian 

The objective of this section is to explain how the interpretation of the Stokes lines used in 
section 2, which led to the rule for the selection of transition sequences, may be justified. 
This requires the use of asymptotic series to reduce the error of the dominant solutions as 
far as possible. It turns out that this error can be reduced below the magnitude of the 
subdominant solution everywhere, except in the vicinity of the Stokes lines. This idea was 
introduced by Stokes [13] in a discussion of the Airy function and amplified in a book by 
Dingle [14]. It has been applied successfully in a variety of forms and to a variety of physical 
problems by Berry [15-17]. To date, all of the applications to differential equations have 
involved problems in which only two equations are coupled together, e.g., two-state problems 
in adiabatic theory, or one-dimensional semiclassical problems, involving two channels (left 
and right propagating waves) . 

The approach used in these papers must be generalized to cover the multi-level problem. 
The reason is that the interpretation of the Stokes lines depends on studying the divergence 
of an asymptotic series. This divergence is usually determined by the branch point 'closest' 
to the real axis, in the sense of having the smallest value of the imaginary part of the action, 
llmS"!. A means must be found to ehminate the divergence associated with branch points 
close to the real axis in order that the divergence associated with other branch points may 
be revealed, allowing their Stokes lines to be interpreted. This is achieved by 'projecting 
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out' a subspace of the Hilbert space spanned by states having branch points close to the 
real axis. 

Various approaches were tried and found to be unsatisfactory. The method wc use is adapted 
from an approach introduced by Berry [18]. Instead of constructing an asymptotic series for 
the wavefunction, the objective is to construct a sequence of 'renormalized' Hamiltonians 
having the same dynamics as the original problem and for which the adiabatic approximation 
is successively more accurate. The k^^ Hamiltonian of this sequence has off-diagonal matrix 
elements of O(e^). If the off-diagonal elements were to approach zero as A; — > oo, the adiabatic 
approximation would be exact. In this case there would be no non-adiabatic transitions. 
This implies that the pref actors of the 0{e'') terms must diverge as /c — > oo. This divergence 
and its consequences are discussed in section 4. In the remainder of this section we explain 
how the sequence of renormalized Hamiltonians is constructed. The approach follows that 
of Berry [18] quite closely, except that a subspace of the spectrum is 'projected out' of the 
renormalization procedure, leaving its elements C)(e°). 

B. Renormalization excluding a subspace 

We consider a case where there is a branch point close to the real axis in a subset of the 
spectrum characterized by a set of state indices {P}- The projection operator P{t) for the 
corresponding sub-space of the Hilbert space is 

p{t)^^\Mt)){Ut)\ ■ (3.1) 

neP 

The complementary set of state indices will be termed {Q} and its corresponding projection 
operator designated Q = I — P. Provided the singularities of the matrix elements of the 
Hamiltonian are sufficiently far from the real axis, the singularities of the projection operator 
P are determined by the branch points where an eigenvalue from {P} becomes degenerate 
with one from {Q}. The projection operator P{t) and the projected Hamiltonian Hp — 
PHP are analytic inside a strip Ep symmetric about the real axis and bounded by these 
conjugate branch points. This is true despite the fact that the eigenstates in (3.1) have 
singularities closer to the real axis due to branch points between states in {P}. 
It will be useful to have available a representation of the projection into the {P} subspace 
explicitly constructed from analytic quantities. To this end a set of states \Xn{t)), n G {P} 
will be constructed having no singularities inside the strip Ep and satisfying 

(Xn(^)IXn'(^)) = Snn' , P{t)\Xn{t)) = \Xn{t)) , lim \Xn{t)) = |0n(±Oo)) . (3.2) 

t—*±oo 

This set of states could be constructed exphcitly in one of several ways. As an example, 
consider the smooth interpolation 

|(|0„(+oo)) + |0„(-oo))) + |(|0„(+oo)) - |0n(-oo))) tanh(et) . (3.3) 

These state vectors could then be projected into the subspace by multiplication by P and 
then orthogonalized by a Gram-Schmidt procedure. 
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Consider the unitary operator 

Mt) = E IXn(t))(Xn(-00)| + ^ \Ut)){U-00)\ (3.4) 
neP neQ 

which generates the states at t from the states eikt — — oo. In order to simplify notation, the 
states |0n(O) '^ill be used to denote the states \Xn{t)) for n e {P} and |0n(^)) for £ {Q} 
so that (3.4) is 

Uo{t)-E\Mt)){k{-oo)\ . (3.5) 

n 

Following Berry [18], we introduce a new representation \ipi{t)) of the solution of the time- 
dependent Schrodinger equation defined by 

m)) = Uo{t)\Mt)) ■ (3-6) 
This renormalized wavefunction satisfies the Schrodinger equation 

ihdtliJi) = Hi\iPi) (3.7) 
with the renormalized Hamiltonian given by 

H,{t) = U^it)H{t)Uoit) - iniJl{t)dMt) . (3.8) 

The matrix elements of H\ are conveniently evaluated in the basis formed by the eigenvectors 
at t — > — oo 

i/W = (0„(-oo)|^i|(/)^(-Oo)) = {4>n{t)\H\Mt)) -H4>n{t)\dtMt)) ■ (3-9) 

Three different cases arise. In the case where n and m are both in {Q} 

^ii^ = <^n.K(t)-m^^^^|J|^^ (n,me{g}) (3.10) 

and for n, m both in {P} 

^^iiil=(Xn(t)|^(t)|Xm(t))-iMXn(t)|atXm(t)) , (n, m G {P}) . (3.11) 
The case n e {P}, m e {Q} requires some discussion. The state \dt(t)m{t)) may be written 

\dt<i>m{t)) = E «n|Xn(t)) + E «n|0n(t)) ■ (3.12) 
ne{P} ne{Q} 

Differentiating the Schrodinger equation {H — £'m)|0m) = with respect to time, then 
multiplying by {xi\ gives the following linear equations for the a„ with n e {P} 

E -^'"^^ = ' (3-13) 

ne{P} 

where 

An = E^Sin - {Xl\PHP\Xn) , = {XlW^m) ■ (3.14) 
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It is useful to define a matrix G that is tfie inverse oi D — {Din\i and a corresponding 
operator G 

G{t)= J2 \xi{t))Gin{t){xn{t)\ = {E^-PHP)-' . (3.15) 

Z,ne{P} 

Using the fact that G = {Gin} is the inverse of D, one finds the solution of (3.13) is 

n€{P} 

= {Xl\{Em-PHPV\Xn){Xn\dtH\(t>m) . (3.16) 
ne{P} 

It follows that 

H^rn = -ihixnimEm ' PHP)-' PdtH\4>m{t)) , (u G {P} , m G {Q}) . (3.17) 

In summary, the matrix dements of the renormalized Hamiltonian are given by the expres- 
sions in (3.10), (3.11), and (3.17). Notice that all of the off-diagonal elements are 0(e), 
except those in the n, m G {P} block, which are 0(1), and that H^^^t) is diagonal as 
t — > ±00, because the original Hamiltonian is constant in both limits. 

C. Iteration of the transformation 

The renormalization transformation defined in the last section can be iterated by writing 
successive wavefunctions and Hamiltonian operators as 

\Mt)) = UkitMk+iit)) , Hk+i = UlHkUk - ihUldA ■ (3.18) 

Note that the transformation of the Hamiltonian is valid for any choice of the unitary 
operator tj^. The operators tJk{t) will be defined by 

Uk{t)-Y.\~€\t)){^^::\-<^)\ (3.19) 

n 

where 

{H,-E^J:Wn^) = Q, (nG{Q}) 

Alxf ) = Ixf ) , (xf (^)lxi'^(^)) = 8nn' , {n,n' G {P}) . (3.20) 

Here the new projection operator P^. is defined in analogy to (3.1) using the new instan- 
taneous eigenstates for n G {P}- The notation |0^f^(t)) is used to denote the states 
IxiiHt)) for n G {P} and |0^f-'(t)) for n G {Q}- Since both iJJj^(±oo) are diagonal, we 
have \(j)!^\±oo)) = |0„(— cxd)). The IXuK^)) are chosen so that they are analytic in the strip 
Epj^^ and close to the |0„(— 00)) as t — > ±00. As in the last section for the k — case, 
these states can be generated by acting on linear combinations of the |0n(±cxD)) with the 
projection operator Pk{t), then using a Gram-Schmidt procedure to create orthonormalized 
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states. It is still most convenient to evaluate matrix elements of the Hamiltonian H^. with re- 
spect to the states oo)). Generalizing the equations in section 3.2, the matrix elements 
H^m^^ = (0„(-oo)|^fc+i(t)|0m(-oo)) are given by 

HLT> - K^E^^Ht) - _ („,„ e {«}) (3.21a) 

= (xi''Ht)I^.WIxL'Ht)) {n,m G {P}) (3.216) 

= -i;i(x^f)(t)|(E« - PkHkPk)-'PkdtHk\<Pl^\t)). [n e {P},me {Q}) (3.21c) 

After k stages of iteration all matrix elements are 0{e''), except the diagonal elements 
and all of the elements of the PP block, which remain 0(1). Each Hamiltonian is an 
exact representation of the dynamics of the original problem, and the transition am- 
plitudes could, in principle, be obtained by integrating the Schrodinger equation with 
any Hamiltonian in the sequence. At each stage the H^^^{±(X)) are diagonal and hence 
|0^'^-'(±oo)) = |0„(— oo)). The amplitude for transition from state n to state m is there- 
fore given by (0„i(oo)|^n(oo)) = (0,„(-oo)|^/'fc,„(oo)), where |V^„(t)) and |^fc,n(t)) are the 
wavefunctions obtained by propagating the initial state |0„(— oo)) with Hamiltonians H 
and Hk, respectively. If the 0{€^) matrix elements vanished as — > oo there would be no 
non-adiabatic transitions. Therefore, the sequence of renormalized Hamiltonians is expected 
to have the typical behavior for terms of an asymptotic series, in that, although the matrix 
elements decrease for small k, at sufficiently large k they diverge because of a faster than 
exponential growth of the prefactors. The value of k for which the largest of the small matrix 
elements has smallest magnitude will be denoted by A;*(e). This divergence will be examined 
in greater detail in section 4. 

For small e and for large values of k which are not too large, the off-diagonal matrix elements 
outside the PP block are very small. Thus, for states in {Q}, the approximations 

~ I0n(-oo)) , Ei'\t) E,,(t) , (e«l, l^k<k*, ne{Q}) (3.22a) 

are vahd, while for states in {P} we are at hberty to choose 

Ixi'H^)) ~ K(-oo)) , (e«l, 1«^<F, ne{P}). (3.226) 

Prom this, along with (3.11) and (3.21b), it follows that for small e 

{xi'\t)\H,{t)\x^^\t)) ~ {xnimmUt)) , {n,m G {P}) . (3.23) 

When these approximations are valid, equations (3.21) simplify considerably. The basis 
states are approximately time independent, so the derivatives with respect to the Hamilto- 
nian become derivatives with respect to matrix elements 

H^nt'^ - KmEn{t) - ^If ^^^^^L (1 ' ^nm) , (n, m G {Q}) (3.24a) 

r^H^lit), (n,mG{P}) (3.246) 

~ -\h G^idAH^^l] , (n G {P}, m G {Q}) . (3.24c) 

ie{P} 
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To summarize: equations (3.21) describe an exact renormalization of the Hamiltonian, which 
for sufficiently large (but not too large) k, has the effect of reducing the magnitude of all 
of the off-diagonal elements outside the PP block. Equations (3.24) are an approximate 
implementation of this iteration valid for sufficiently small e. 

IV. INTERPRETATION OF THE ASYMPTOTIC SERIES 

All of the renormalized Hamiltonians Hk describe the same dynamics. In particular 
|(0m(— oo)|-?/'fc „(cxd))P, where \ipk,n(t)) is the wavefunction propagated under the Hamilto- 
nian Hk from the initial state |0„(— oo)), is the transition probability from state n to state 
m and is independent of k. Now the viewpoint of section 2, which associates transitions 
between pairs of levels with Stokes lines, will be confirmed here by showing that when k 
is suitably large, the Hamiltonian Hk{t) on the real axis is greatest at crossings of Stokes 
lines [17]. Section 4.1 examines the case where there is no projected subspace: a direct 
generalization of the result for two level systems. Section 4.2 considers the extension to the 
case where there is a projected subspace, showing that the Stokes lines can still have the 
same significance for branch points not close to the real axis. Section 4.3 describes how 
these results can be used to calculate a transition probability involving a sequence of branch 
points. 

A. Form of the Hamiltonian at large order, excluding projected subspace 

This section discusses the form of the matrix elements of the renormalized Hamiltonian in 
the limit where e <S 1 and 1 <^ k < k*. Equation (3.24a) which determines the evolution of 
off-diagonal matrix elements in the QQ block, is an iteration of the form 

Hk+i(t) = -^d,H,{t) (4.1) 

where Hi{t) and A(t) are given functions. From observations on asymptotic series discussed 
by Dingle [14] it is well known that for large k, solutions of (4.1) may be obtained in the 
form 

H,{t) ~ (-l)^Ar(A; + 7)[F(i)]-('+^) (4.2) 

where A and 7 are constants, and r{x) is the Gamma function satisfying T{x + 1) = xT^x) = 
x\. Substitution of (4.2) into (4.1) shows that the dimensionless scalar function F{t) is 

F{t) ^ y^Jt' A{t') , (4.3) 

and the lower limit t* will sometimes be shown explicitly as F{t,t*). The possible values 
of t* are determined by realizing that an acceptable solution of (4.1) should be regular at 
any point for which A~^(i) and Hi{t) are regular. The solution (4.2) is clearly divergent 
for k + J > when t ^ t*. This observation implies that t* must correspond to a singular 
point, where at least one of A-\t) and Hi{t) has non-analytic behavior. In the context of 
adiabatic theory t* corresponds to a branch point singularity. 
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In general, the solution at time t is dominated by the approximate solution (4.2) attached 
to a singular point f* for which the modulus of F(t,t*) is smallest. Since equation (4.1) is 
linear, a superposition of co-dominant solutions of the form (4.2) may be necessary. In the 
case where the functions Hi{t) and A~^(t) are analytic and real valued on the real axis, the 
singular points occur as complex conjugate pairs. In this situation the solution on the real 
axis is a superposition of solutions with conjugate singularities. 

In our application where A = {Em — En), the real part of F is constant along the real axis. 
Thus, on the real axis, the magnitude of the solution (4.2) is largest at the point where the 
Stokes line, defined by lmF{t,t*) = 0, crosses. For large k the magnitude decays rapidly, 
with an approximately gaussian form on either side of this point. When the two conjugate 
solutions are combined, the result is a real function with oscillations within an approximately 
gaussian envelope. 

The constants A and 7 can be determined - although they are not required for appli- 
cation of the theory. Both constants are obtained by considering the behavior of the 
functions in the neighborhood of the singularity. In the application of these results 
to adiabatic theory, equation (3.10) requires that the function Hi{t) plays the role of 
{(l)n{t)\dtH\(t)mit)) / (Emit) — En{t)) , which scalcs as c(t — t*)~^ in the vicinity of a branch 
point singularity at t*. In the vicinity of the branch point, the form of the successive Hf^{t) 
is dominated by the components of Hi{t) and A~^(i) which are most strongly divergent as 
t^t*. It is therefore appropriate to take Hi{t) ^c{t-t*)-^ and A(t) = Q;(i - r)^^ ^his 
leads by direct application of (4.1) to a general form for the Hk{t) vahd in the neighborhood 
of t* 

H,{t)^C,{t-tT^''-'^^\ C,^i^-^^C,. (4.4) 

za 

The iteration for the coefficients Ck has a solution 

Comparison of this result with (4.2) shows that they are consistent in the limit k ^ 00, and 
that 7 = -i and A = -(f )ic/r(|). 

Applying these results to the renormalization of the QQ block of the Hamiltonian (equation 
(3.24a)), shows that for large k the off-diagonal matrix elements are largest for the element 
H^l!^l for which \F\ is smallest. On the real axis F = iSnm/T^, so this corresponds to the matrix 
element with the least |Im5'„m|- Furthermore, each of these matrix elements is largest in 
the vicinity of the point where the Stokes line intersects the real axis. The most significant 
process in the dynamics of is therefore a transition between levels n and m localized in 
time at the point where the Stokes line crosses the real axis. The use of the renormalization 
scheme therefore confirms the interpretation that the transition occurs at the Stokes line. 
The results are consistent with the conclusions of section 2, and we will assume that the 
transition probability is given by (2.14) and (2.15). 

B. Form of the Hamiltonian at large order, including projected subspace 

In the case where there is a sub-space projected out in the manner discussed in section 3, 
the equation for the iteration of the off-diagonal matrix elements in the QQ block is exactly 
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of tfie form considered in section 4.1 above. The case of matrix elements in the PQ and QP 
blocks requires a different treatment. The equation for iteration of the these elements of the 
Hamiltonian is of the form (3.24c), which we write symbolically as 

Hk+i{t)^-ihG{t)dtHk{t) (4.6) 

where Hk is a column vector of dimension A/p, with elements H^J^, and G = {Gni} = 
{{Xn\{Em — PHP)^^\xi)}, a matrix of dimension A/p x A/p. Note that the index m is fixed 
and not shown in (4.6). Following (4.2), it might be anticipated that a solution of (4.6) is 
of the form 

H, = (-l)^r(A: + 7)[M(t)]-^'+^^i , (4.7) 

where M{t) is a matrix of dimension A/p x A/p and A a constant column vector of dimension 
A/p. We proceed to show that (4.7) is the solution to (4.6) to leading order in k~^. Let 
M{t) = X^^{t)F{t)X{t), where F{t) is diagonal. Then the time derivative may be simplified 
using the following approximation valid to leading order in k~^: 

dt[M{t)r'' = -kX-\t)dtF{t)[F{t)]-^''^^^X{t) + 0{1) . (4.8) 

Using this approximation one finds that (4.7) is the solution to (4.6) to leading order in k~^ 
so long as —ihXGX'^dtF is the identity matrix. It follows that XGX~^ must be diagonal 
since F is assumed to be diagonal. The matrix X{t) is therefore the similarity transformation 
that diagonalizes G{t), and the matrix F(t) and its elements Fn{t) are 

F(t) = ^ r dt' ^{t') G(t')-' X-\t') 

Fn{t) = ^ £ dt' [Emit') - Er^it')] . (4.9) 

Substitution into (4.7) now gives the required approximation. These results are a natural 
generalization of the results for the simpler case described by the iteration (4.1). A similar 
argument to that given in section 4.1 forces t* to be the position of a singularity, where a 
pair of eigenvalues, one each from the P and Q subspaces, become degenerate. 
It is worth noting that these expressions can be simplified by making use of freedom available 
in the specification of the states \xn{t))- The states \Xn{t)) may be chosen to coincide with 
the eigenstates |0„(t)) at isolated positions in the complex time plane, so that the similarity 
transformation X(t) is the identity at that point. Specifically, these states may be chosen 
to coincide at the point where a Stokes line crosses the real axis. The comments in section 
4.1, indicating how the renormalized Hamiltonian is concentrated on Stokes lines, can then 
be apphed equally well to this case: from (4.7) and (4.9) it can be seen that the dominant 
term in the high-order renormalized Hamiltonian comes from the branch point with the 
smallest value of |lmS'„m|) where at least one of n, m is in the Q subspace. If X{t) is 
the identity where the corresponding Stokes line crosses the real axis, the largest matrix 
elements determine the transition between states n and m. 

C. Interpretation of the renormalized Hamiltonian 

As an example of how these considerations can be used to demonstrate the validity of the 
topological rule, consider their application to the case of two successive transitions Uq — > rii 
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followed by rii —^n2. It will be assumed that the (rio, ^i) branch point is closest to the real 
axis and that the (ni, 77.2) branch point is next closest. The subspace P is then chosen to be 
that spanned by levels no, rii. The states \Xnoii)) and |Xni(^)) are chosen to correspond to 
the eigenstates \(t>noif)) and at the point ts where the Stokes line from the (ni,n2) 

branch point crosses the real axis. When the order k is sufficiently large, the matrix elements 
that determine transitions from the P subspace are largest for the transition between levels 
rii and n2, and are concentrated at ts- The general arguments given in section 2 indicate 
that this transition probability is P„j^„2 ^ exp[— 2 Im S'„j ^j/^]. 

The probability for making a transition from state Uq to n2 via the intermediate state rii 
therefore depends on the dynamics within the PP sub- block at times earlier than ts- The 
matrix elements within this sub-block have been left approximately unchanged by iteration 
of the renormalization procedure. Any suitable procedure can be used to calculate the 
probability for the transition no fii occurring before the time ts- In the limit e — > 0, the 
most convenient procedure is of course to use adiabatic theory that predicts the transition 
will occur on the Stokes line attached to the (no,ni) branch point with the probability 
Pno^ni ~ exp[— 2 Im ^no.nj//?.]. Thc ovcrall transition probability for the no n2 transition 
going through the intermediate state rii is then P„(,^„2 = P„(,^„^P„^^„2, if the branch point 
(ni, n2) lies to the right of the Stokes line from the (no, ni) branch point. In the other case, 
it is determined by other branch points. 



V. NUMERICAL ILLUSTRATION 

The scenario described above was tested numerically on a model Hamiltonian H{X{t)) of 
the following form: 

^ = cos(X(r))i7i + sin(X(r))^2 (5.1) 

A:(r) = atanh(r/a) (5.2) 

where r = et, and Hi and H2 are two independent random square matrices of dimension 
A/", drawn from the gaussian orthogonal ensemble (GOE). The GOE ensemble consists of 
real, symmetric matrices with independently gaussian-distributed elements with mean and 
variance given by [20,21] 

= 0, (//J) = (1 + 5,,.) . (5.3) 

The choice of Hamiltonian is arbitrary, so long as its spectrum is nondegenerate for all real 
r, and its eigenvalues and eigenstates become time independent asymptotically as r — > ±oo- 
The function X{t) was chosen to fulfill the latter requirement. The model (5.1) was used 
because, for J\f ^ 1, its spectral properties are representative of those of generic time- 
reversal invariant multi-level systems [22]. The distance of the branch point singularities 
from the real axis scales as M~^^'^, whereas the singularities of the matrix elements in the 
complex r plane have a distribution of positions independent of the matrix dimension M - 
It follows that in typical physical applications branch points will lie closer to the real axis 
than other singularities. 

Numerical calculations were performed to determine the 'exact' or 'empirical' transition 
probabilities Pn^m for 100 Hamiltonians of the form (5.1), all with dimension J\f — 6. In each 
case all J\f{J\f — 1) = 30 transition probabilities were computed. The scale factor a = 2 was 
assumed in all calculations. The probabilities were obtained by numerical time integration 
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of the Schrodingcr equation, using a standard fourth order Runge-Kiitta algorithm with 
arithmetic accurate to 14 decimal places. Rather than comparing the transition probabilities 
themselves, we compared (the imaginary part of) the actions times the adiabatic parameter 
e in units of Planck's constant: 

I Im Snm I 1 

The integration was performed between = —25 and Tf = +25. We checked that the 
transition probabilities are insensitive to further increasing the range of r. These calculations 
were done for a large set of adiabatic parameters between e = 0.001 and e = 0.5. We observed 
that in most cases the actions converged to a limit as e — > 0. In other cases, we estimated 
the limit by extrapolation with a polynomial. In many cases numerical roundoff error would 
make the results unreliable for small e. In these cases we assumed the best value for the 
action corresponded to the smallest value of e for which roundoff errors were not significant. 
Typically, roundoff errors were significant when the calculated result obeyed Pn-^m < 10~^^. 
The adiabatic calculations were performed in the following way. First, the approximate 
locations of the branch points t*j were determined by a search for near- degeneracies on 
a grid in the complex r plane. The rectangular region bounded by —6 < Re r < 6 and 
< Im r < 2 was usually found to include all the relevant branch points. Next, the locations 
of these branch point candidates were refined by a version of the Newt on- Rap hson method 
adapted to finding square root branch points of the form — Ej oc — rlj. The Stokes 
and anti-Stokes lines were plotted for several realizations of the random Hamiltonian, and the 
allowed transition sequences were determined. The Stokes and anti-Stokes lines attached to 
the branch points were traced by evaluating a sequence of short steps 6t along their lengths. 
If the Stokes line emerging from a branch point between levels i and j was found to pass 
through the point t^, the next point was obtained from r^+i — Tk + Srk where from (2.6) 
^Tfe = ei/[Ei(Tk) — Ej(rk)], and £ is a small real number. Increments of the anti-Stokes lines 
were determined by an analogous approach. Finally, the 'adiabatic' or 'theoretical' branch 
point actions were computed numerically according to 



Im dr {E,{t) - E,{t)) . (5.5) 

•^(Rer*.,0) 



The theoretical actions were in excellent agreement with those determined empirically from 
(5.4). We found that in every case the allowed transition sequence could have been deter- 
mined by a simple empirical rule, namely, that the real parts of the branch points should be 
in ascending order. Accordingly, the branch point data were ordered with respect to ReT*j. 
From this ordered table of branch point locations and branch point actions, all possible 
transition sequences from the initial state n to the final state m, in which the state indices 
increased (or decreased) monotonically, were considered. According to the rule given in 
section 2.5, the overall action A„m was taken to be the least sum of (the imaginary part of) 
the branch point actions over these possible transition sequences. These calculations of the 
overall 'theoretical' actions were automated. 

Examples of the calculations are presented in the tables and figures. Figures 6 and 7 show 
plots of the energy eigenvalues for real time for two examples of the 100 sample Hamiltonians. 
The branch point data corresponding to these two cases are listed in tables 1 and 2. Observe 
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that the avoided crossings in these figures correspond to small values of the branch point 
actions in the tables. Figures 8 and 9 show the Stokes lines that cross the real axis for all 
branch points involved in some transition sequence. The 'theoretical' transition sequences 
and actions derived from these data are shown in tables 3 and 4. The 'empirical' results 
obtained using (5.4) are shown there for comparison. Note that the fractional difference 
between the empirical and theoretical results based on our adiabatic theory is typically 
around 1% or smaller, although occasionally as high as 2%. The average fractional difference 
between the empirical and theoretical actions over all 3,000 data points was ~ 1.6%. The 
data are consistent with the hypothesis that the transition probabihties are given by (1.1) 
with C — 1, and the action Snm given by the topological rule of section 2.5. 

VI. CONCLUDING REMARKS 

This paper makes three novel contributions to understanding the behavior of the transition 
probabilities for multi-level systems. First, we have suggested a general rule for determin- 
ing the combinations of branch points that give allowed transition sequences, based upon 
an assumption that the transitions occur when a Stokes line, whose dominant solution is 
occupied, is crossed. Second, we have indicated a general approach to interpreting all of the 
Stokes lines in a multi-level system, by using a projection technique to selectively remove 
the most divergent contributions from the renormalized Hamiltonian. And finally, we have 
verified the rule by means of extensive numerical experiments. 

The results of sections 3 and 4 support the rule suggested in section 2, but they do not 
constitute a proof. Further work must be done to establish how the high order renormal- 
ized Hamiltonian can be used to calculate transition probabilities directly in the case of 
multi-level systems. This might involve studies of how the transition probabilities could be 
determined by applying time-dependent perturbation theory to the high order renormalized 
Hamiltonians as t increases along the real axis. This has been successfully applied by Berry 
[17] for the case of two level systems. An alternative approach would be to investigate an- 
alytic continuations of the adiabatic solutions of the renormalized Hamiltonian away from 
the real far as the branch points. 

There is a theoretical difficulty that must be resolved concerning the interpretation of the 
Stokes lines. According to the interpretation discussed in section 4, a transition from level 
no to rii occurs on crossing a Stokes line S{no,ni), where the matrix elements H^J,^^{t) 
are greatest for real t. A subsequent transition from state Ui to n2 could then occur if 
the Stokes line S{ni,n2) were crossed at a later time. Presumably the allowed transition 
sequence would then be determined by the order in which the Stokes lines cross the real 
axis. However, the path in the complex time plane along which the Schrodinger equation 
is integrated can be deformed away from the real t axis. A problem may arise if the two 
Stokes lines cross: a transition allowed along one path would then be forbidden along an 
equally valid path. A rule based on the order in which Stokes lines cross the real axis may 
therefore predict different transition sequences from the rule in section 2. We have not yet 
found a totally satisfactory resolution of this problem. We can, however, remark that in our 
investigations we found no examples where the predictions of the two possible rules were 
different. It is not clear whether it is impossible to find an example in which the predictions 
differ, or whether such cases occur with very low probability within our ensemble. 
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FIGURES 



Figure 1. The exponent characterising a non-adiabatic transition in a two-level system is 
obtained by integrating the energy along a path enclosing a branch point in the complex t 
plane where the levels become degenerate. 

Figure 2. (a) Energy levels as a function of real time for a three level system. The regions 
where the curves approach each other closely are called 'avoided crossings', (b) The avoided 
crossings correspond to branch points close to the real axis. 

Figure 3. (a) The sign of an eigenvector is reversed upon making two circuits about a branch 
point, (b) Illustrates this change of sign for the Landau-Zener Hamiltonian by considering 
the eigenvector transported smoothly around a circuit taken to infinity in the upper-half 
plane. 

Figure 4. Stokes lines (Sn), anti-Stokes lines (An), and branch cut (BC) associated with a 
single branch point. 

Figure 5. Stokes and anti-Stokes lines associated with a pair of branch points between 
different pairs of levels. The boundary condition is that only level uq is occupied as t — > — oo. 
Case (a) allows a transition to level n2 stst ^ oo. Case (b) does not allow such a transition. 

Figure 6. Energy levels for example I of the model Hamiltonian introduced in section 5. 

Figure 7. Energy levels for example II of the model Hamiltonian introduced in section 5. 

Figure 8. Branch points and Stokes lines intersecting the real axis for example I. The pairs 
of integers indicate which energy levels become degenerate at each of the branch points. 

Figure 9. Branch points and Stokes lines intersecting the real axis for example II. 
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TABLES 

Table 1. Branch point data for example I. 
Table 2. Branch point data for example II. 

Table 3. Theoretical transition sequences and actions compared with empirical values for 
example 1. 

Tabic 4. Theoretical transition sequences and actions compared with empirical values for 
example 11. 
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branch 


State 


Re Tij 




point label 


indices i,j 


a 


5,4 


-3.051296 


2.637619 


b 


4,3 


-2.739876 


.181812 


c 


2,1 


-1.971837 


2.918432 


d 


3,2 


-1.791580 


.049935 


e 


6,5 


-1.446037 


.101999 


f 


4,3 


-1.195295 


.119999 


g 


5,4 


-.729159 


.031601 


h 


2,1 


-.610604 


.024514 


i 


6,5 


-.440841 


.267925 


i 

J 


4,3 


-.223932 


.256871 


k 


3,2 


.296616 


.071114 


1 


5,4 


.668725 


.122366 


m 


5,3 


.678210 


.480203 


n 


2,1 


.790463 


.373673 


o 


4,3 


.884734 


.166006 


P 


3,2 


1.242660 


.304654 


q 


6,5 


1.609035 


2.043553 


r 


4,3 


1.647808 


.072687 


s 


5,4 


2.528592 


.090388 


t 


2,1 


3.066494 


.452725 



table 



Im T 
















i I' 

— > 




E=E, 


^=^2 Re T 



fig.1 



branch 
point label 


state 
indices i,j 


ReXii 


''J 


a 


4,3 


-2.815384 


.243081 


b 


6,5 


-2.338206 


1.807943 


c 


2,1 


-2.134465 


1.915952 


d 


3,2 


-2.065168 


.482255 


e 


6,4 


-2.030557 


5.220234 


f 


4,3 


-1.809887 


.327162 


2 
o 


2,1 


-1.745756 


.514480 


h 


5,3 


-1.655970 


3.443324 


i 


3,2 


-1.291675 


.792717 


i 

J 


5,4 


-1.236497 


.019950 


k 


4,3 


-.564900 


.028096 


1 


6,5 


.096006 


1.027531 


m 


5,4 


.351965 


2.323682 


n 


2,1 


.904734 


.936584 


o 


3,2 


.965953 


.027680 


p 


3,1 


.982541 


1.177381 


q 


4,2 


1.037342 


1.035606 


r 


6,5 


1.086503 


.476258 


s 


2,1 


1.373743 


.645514 


t 


4,3 


1.579084 


.085708 


u 


5,4 


1.986832 


.465658 


V 


6,3 


2.181664 


3.069844 


w 


4,3 


2.191573 


.428753 


X 


6,5 


2.227829 


.770383 



table 



Initial 


Final 


Transition 


Km 


Km 


Difference 


Fractional 


state n 


state m 


sequence 


adiabatic 


empirical 




difference 


1 


2 


h 


.024514 


.024972 


-.00046 


-.01833 


1 


3 


hk 


.095628 


.097049 


-.00142 


-.01463 


1 


4 


hkr 


.168315 


.170316 


-.00200 


-.01174 


1 


5 


h k r s 


.258703 


.261363 


-.00266 


-.01018 


1 


6 


h k m q 


2.619384 


2.630638 


-.01125 


-.00428 


2 


1 


h 


.024514 


.024946 


-.00043 


-.01732 


2 


3 


d 


.049935 


.050603 


-.00067 


-.01319 


2 


4 


dr 


.122622 


.123828 


-.00121 


-.00974 


2 


5 


dfg 


.201535 


.204990 


-.00345 


-.01685 


2 


6 


dfgi 


.469460 


.473734 


-.00427 


-.00902 


3 


1 


dh 


.074449 


.075558 


-.00111 


-.01467 


3 


2 


d 


.049935 


.050611 


-.00068 


-.01336 


3 


4 


r 


.072687 


.073213 


-.00053 


-.00719 


3 


5 


fg 


.151600 


.152904 


-.00130 


-.00853 


3 


6 




.419525 


.423847 


-.00432 


-.01020 


4 


1 


bdh 


.256262 


.257967 


-.00171 


-.00661 


4 


2 


fk 


.191113 


.193036 


-.00192 


-.00996 


4 


3 


r 


.072687 


.073212 


-.00053 


-.00717 


4 


5 


2 
o 


.031601 


.032305 


-.00070 


-.02181 


4 


6 


g i 


.299526 


.302137 


-.00261 


-.00864 


5 


1 


g i k n 


.733259 


.736357 


-.00310 


-.00421 


5 


2 


g i k 

O J 


.359586 


.362679 


-.00309 


-.00853 


5 


3 


g r 


.104288 


.105613 


-.00133 


-.01255 


5 


4 


g 


.031601 


.032305 


-.00070 


-.02179 


5 


6 


e 


.101999 


.102708 


-.00071 


-.00691 


6 


1 


egjkn 


.835258 


.846646 


-.01139 


-.01345 


6 


2 


egjk 


.461585 


.465363 


-.00378 


-.00812 


6 


3 


egr 


.206286 


.208215 


-.00193 


-.00926 


6 


4 


eg 


.133599 


.134977 


-.00138 


-.01021 


6 


5 


e 


.101999 


.102680 


-.00068 


-.00664 



table 3 



Initial 
state n 


Final 
state m 


Transition 
sequence 


Km 

adiabatic 


Km 

empirical 


Difference 


Fractional 
difference 


1 


2 


2 
o 


.514481 


.516859 


-.00238 


-.00460 


1 


3 


o 


.542161 


.544292 


-.00213 


-.00392 


1 


4 


got 


.627869 


.630412 


-.00254 


-.00403 


1 


5 


g o t u 


1.093527 


1.098524 


-.00500 


-.00455 


1 


6 


g o t U X 


1.863910 


1.870412 


-.00650 


-.00348 


2 


1 


g 


.514481 


.516342 


-.00186 


-.00360 


2 


3 


o 


.027680 


.028112 


-.00043 


-.01535 


2 


4 


ot 


.113388 


.114395 


-.00101 


-.00880 


2 


5 


o t u 


.579047 


.581847 


-.00280 


-.00481 


2 


6 


dfjr 


1.305625 


1.310825 


-.00520 


-.00397 


3 


1 


o s 


.673194 


.677931 


-.00474 


-.00699 


3 


2 


o 


.027680 


.028111 


-.00043 


-.01532 


3 


4 


k 


.028096 


.028485 


-.00039 


-.01365 


3 


5 




.263031 


.263790 


-.00076 


-.00288 


3 


6 


a i r 


.739289 


.742328 


-.00304 


-.00409 


4 


1 


k o s 


.701290 


.708465 


-.00717 


-.01013 


4 


2 


ko 


.055776 


.056596 


-.00082 


-.01448 


4 


3 


k 


.028096 


.028485 


-.00039 


-.01365 


4 


5 


i 

J 


.019950 


.020313 


-.00036 


-.01786 


4 


6 


i r 

J 


.496208 


.498257 


-.00205 


-.00411 


5 


1 


i ko s 


.721240 


.724360 


-.00312 


-.00431 


5 


2 


i k o 

J 


.075726 


.076910 


-.00118 


-.01539 


5 


3 


jk 

J 


.048046 


.048794 


-.00075 


-.01534 


5 


4 


i 

J 


.019950 


.020312 


-.00036 


-.01784 


5 


6 


r 


.476258 


.479487 


-.00323 


-.00673 


6 


1 


b j k o s 


2.529184 


2.532333 


-.00315 


-.00124 


6 


2 


bjko 


1.883670 


1.888210 


-.00454 


-.00240 


6 


3 


r u w 


1.370669 


1.376297 


-.00563 


-.00409 


6 


4 


r u 


.941917 


.944955 


-.00304 


-.00322 


6 


5 


r 


.476258 


.478059 


-.00180 


-.00377 



table 4 
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